COMPUTING DULAC'S MAPS OF (ALMOST EVERY) 
NON-DEGENERATE SINGULARITIES 



LOIC TEYSSIER 



Abstract. We study the complex Dulac map for a holomorphic foliation near a 
non-degenerate singularity (both eigenvalues of the linearization are nonzero). We 
describe the order of magnitude of the first two terms of the asymptotic expansion 
and show how to compute explicitly those terms using characteristics supported in the 
leaves of the linearized foliation. We perform similarly the study of the Dulac time 
spent around the singularity. These results are formulated in a unified framework 
taking no heed to the usual dynamical discrimination (i.e. no matter whether the 
singularity is formally orbitally linearizable or not and regardless of the arithmetics 
of the eigenvalues) provided the foliation has enough (i.e. two) separatrices. 



1. Introduction 

We consider a germ of a holomorphic vector field at the origin of C 2 

A(x,y)^+B(x,y)^ 

admitting an isolated, non-degenerate singularity at (0,0). In other words the origin is 
the only local zero of the vector field, and its linear part [VA , VB] at this point is a 2 x 2 
matrix with two nonzero eigenvalues, of ratio A G C^o- We make the further assumption 
that the vector field admits the following decomposition in a convenient local analytic 
chart: 

(X) X R = \x— + (1 + R) y— , R G x a C {x, y} 

for some positive integer a satisfying the relation 
(a) Sftfa+-M>0. 

This setting encompasses almost all non-degenerate singularities, including every kind of 
saddle singularities (A < 0), as will be discussed at the end of the introduction. 

Our study is carried out on a fixed polydisc U = pD x rO small enough for the relation 
(R) sup|i?(^)|<l 

to hold. We write Tr the foliation of U whose leaves are defined by the integral curves of 
Xr. This foliation admits two special leaves (called separatrices) each of whose adherence 
corresponds to a branch of {xy = 0}. We write U := U\ {xy = 0}. Outside {x = 0} the 
foliation is transverse everywhere to the fibers of the fibration 

n : (x,y) i — > x. 

Being given (x*,y*) ^ U it is thus possible (under suitable assumptions that will be 
detailed later on) to lift through II a path 7 linking x to 2;* in the foliation, starting 
from the point (re,?/*). The arrival end-point of the lifted path defines uniquely a point 
{x*, y x ) £ II -1 (2:*). This construction yields a locally analytic map from a sub-domain of 
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the transverse disc {y = y*} into the transverse {x = x*}, which is known as the Dulac 
map 

V R : x ^ i — > y x 

of Xr associated to (x*,y*). 
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Figure 1.1. The Dulac map. 

This map is in general multivalued, and its monodromy is generated by the holonomy 
of Tr computed on {x = a;*} by winding around {x = 0}. The Dulac map governs part 
of the dynamics of Tr and has been submitted to an intense study at least in the two 
settings we describe now. 

• The proof of the Dulac conjecture regarding finiteness of the number of limit cycles 
for analytic vector fields in the real plane. Dulac maps are natural ingredients 
of the cross-section first-return map along a singular poly-cycle, whose attractive 
fixed-points correspond to limit cycles. As was noticed by Y.Il'yashenko [I1'85] 
the original «proof» of H. Dulac [Dul23] crucially depends on a lemma which 
turned out to be false. Many powerful, if intricate, tools have been developed in 
the 1980's decade which finally led to a complete proof of Dulac's conjecture. Two 
parallel approaches evolved at the time to analyze the asymptotic expansion of 
the Dulac map: J. ECALLE [Eca92] studied it first formally using trans-series then 
through resurgent summation techniques, while Y. Il'yashenko [IF91] devised 
an argument based on super-accurate asymptotic series. Of course various authors 
contributed to the tale of Dulac's conjecture but the aim of this article is not to 
offer a comprehensive list. 

• The study of the dynamics of a singular foliation through Seidenberg's reduction 1 
process: Dulac maps appear naturally as «corner maps» encoding the transition 
between different components of the exceptional divisor. They measure how the 
different components of the projective holonomy pseudo-group mix together (see 
e.g. [LorlO] for a comprehensive presentation). In that context D. Marin and 
J.-F. Mattei [MM08] proved recently that under suitable (generic) hypothesis 
a germ of a singular foliation is locally incompressible: there exists an adapted 
base of neighborhoods of (0, 0) in which the (non-trivial) cycles lying in the leaves 
of the restricted foliation must wind around the complement of the separatrix 
locus, in the trail of Milnor's theorem regarding holomorphic fibrations outside 

1 According to Seidenberg's algorithm (see [Sei68]) any isolated singularity of a germ of a holomorphic 
foliation T can be «reduced» through a proper, rational map 7r : M — > (C 2 , 0) , where A4 is a conformal 
neighborhood of a tree E := 7r -1 (0,0) of normally-crossing, conformal divisors Pi (C). The pulled- 
back foliation it*T has only isolated, reduced singularities (located on E): either non-degenerate or of 
saddle- node type (exactly one nonzero eigenvalue). 
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the singular fibers. This study is the first step towards a complete analytical and 
topological classification of (generic) singular germs of a foliation. One of the 
main ingredients of their proof is to control the «roughness» of the corner maps 
and elements of the projective holonomy. This roughness can be read in the first 
two terms of the asymptotic expansion of the Dulac map. 

The topic addressed here may appear worthless in view of the formidable work carried 
out by both J. Ecalle and Y. Il'yashenko. Yet it must be noticed that they were 
essentially busy with the real context and, although their techniques are naturally part of 
complex analysis, they focused on small neighborhoods of real trajectories. The originality 
of the present work lies in the study of the analytic continuation of the Dulac map on 
large sectors around the real directions. More precisely our results cover the following 
items, for saddle-like singularities (that is, 5ft (A) < 0): 

• we give a geometrical description of the maximal domain of definition (it contains 
germs of a sector around {x = 0} of arbitrary aperture), 

• we express the Dulac map as an integral (more precisely, as the characteristics of 
the differential 1-form Rd (logs;" 1 /*) along the leaves of Fr), 

• we derive explicit bounds on the remainder of the asymptotic expansion in the 



In particular it should be noted that although the Dulac map's asymptotic expansion can 
be expressed formally 



the expansion converges if, and only if, the foliation is analytically normalizable 2 , as was 
proved by A. MOURTADA and R. MOUSSU [MM97, Proposition 1]. 

We also mention that the result exposed here allows to get rid of a non-necessary tech- 
nical hypothesis in Marin-Mattei's theorem, namely discarding «bad 3 » irrational ratios 
in Seidenberg's reduction of the singularity. Proposition 1.5, stated further down, implies 
that the hypothesis of Lemma 4.4.3 of [MM08] are fulfilled, which was the only depen- 
dency in the chain of proof that needed to be addressed. The other restriction in their 
theorem, namely that no degenerate singularity occurs in the reduction process, is not 
technical and although it can be marginally weakened, must be kept (see [MT13] for ex- 
amples of compressible foliations). Marin-Mattei's theorem now holds on the complement 
of a proper analytic set of the space of singular holomorphic foliations, the one avoiding 
degenerate singularities in the singularity's reduction. 

1.1. Statement of the main results. The Dulac map is actually well-defined on the 
foliated universal covering of the restriction of Fr to U, a complex surface isomorphic 
to the groupoid of paths tangent to Fr with starting endpoint in {y = , x ^ 0}. Our 
first result regards the structure of this space, which turns out to be embedded in the 
universal covering of U. This is another manner of saying that Tr is incompressible in 
U: every non-trivial cycle of a leaf of Tr is a non-trivial cycle of U. 



Situation which arises not so often in the most interesting (quasi-) resonant cases. 
3 A «bad» irrational is unusually well approached by rational numbers. This only happens on a set of 
null Lebesgue-measure, but does nonetheless. Brjuno's explicit arithmetic condition [Brj71], expressed in 
terms of the convergents of A, governs this behavior. 



sectors. 




where 
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Theorem 1.1. Let £ : (z,w) i — > (expz,expu>) be the universal covering oflA. Assume 
that conditions (X) and (R) are fulfilled. Then the foliation £* JFr is regular and each one 
of its leaves is simply- connected. 

We understand now T>r as a function z i— > w z , after having fixed once and for all a 
preimage (z*,w*) G £~ 1 (a;*,j/*). The domain of definition of T>r corresponds to those 
(z, w„) giving birth to a path 7^ (z) tangent to £*Tr with landing endpoint (z*, w z ). We 
state now our main result. 

Theorem 1.2. ylsswme that conditions (X) and (R) hold. 

(1) 2?r is holomorphic on a domain f2» 3 z» satisfying the following additional prop- 
erties for any N € N>o- 

• 7/ 5ft (A) < t/ien f2* is simply- connected and there exists r > r' > suc/i 
£/ia£, for any w* in {5ft (w) < lnr'} and 5ft (z*) < In p ; t/ie domain fi, contains 
at least the line segment 

{»(«) = ln/9, |3?(z-z*)| <7r^V} 

in its adherence. Besides inf 5ft(f2* fl {|3 (z — z*)| < -ZV7r}) > —00. 

• 7/ 5ft (A) < £/iere exists < p' < p such that for every 5ft (z„) < hip' and 
N G N>o i/ie domain f2* contains some infinite half-band 

{5ft (z) < ft', |3(z- Z,)| < tt^} 

iyii/j k' < 5ft (z*). 

(2) .For every z £ £7* we /iai>e 

2? fl (z) = + - + - / Ro£ dz. 

A A AhM 

7n particular T>r (z*) = w„. 

(3) If 5ft (A) < and the condition (a) is satisfied then one has the asymptotic expan- 
sion 

j Ro£ dz = / R o £ dz + o (|z exp ~ z /\\) 

when 5ft (z) tends to —00 with a bounded imaginary part. 

Remark 1.3. In (3) it only makes sense to consider that «5ft (z) tends to —00 with a 
bounded imaginary part» when 5ft (A) < 0, as underlined in (1). This is why we make the 
hypothesis on the real part of A to precise the asymptotic behavior of the Dulac map. 
Although it is a technical artifact that can be overcome, by considering slanted bands 
instead of «horizontal» bands {|9 (z — z*) < ttN}, I surmise that argz must be somehow 
controlled in order to be able to obtain a useful asymptotic expansion as 5ft (z) — > — 00. 

It is possible to carry out the computation of the characteristics associated to the model 
Fq. The exact value of f^ Q r z \ {x n y m ) o £ dz does not offer an insightful interest as such. 
We can nonetheless deduce from it the dominant part of the characteristics f , . Ro£dz, 

which splits into two components: the regular part J * R (exp it, 0) dw, which induces a 
holomorphic function in IA since R (0, 0) = 0, and the resonant part obtained by selecting 
in R only well-chosen monomials. 

Definition 1.4. The resonant support Res (a, A) associated to (a, A) is 

• the empty set if A ^ K<o, 

• otherwise the subset of N 2 defined by 

Res (a, A) := < (n, m) £N 2 : m > , n > a , \nX + ml < — 

2n 
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For G (x, y) = Y, n >o, m>o G n, m x n y m G C {x, y} we denote by G the regular part of G 

Go (a?) := G(x,0) 

and by G Res its resonant part 

G Res (x,y) := Gn. m x n y m . 

(n,m) ERes(a, A) 

It turns out that this support consists indeed of resonant or quasi-resonant monomials, 
according to the rationality of A (Lemma 3.4), which carry without surprise the major 
part of the non-regular characteristics. 

Proposition 1.5. Assume that -ft (A) < and condition (a) holds. Then for any G G 
C {x, y} we have 

( G o Edz = G (0) (z* -z) + (exp z) 
and if moreover G G x a< C {x, y} 

G Res o£dz = O (\zexp- z /x\) 



'7o(*) 

/ (G-G Res -G Q )o£dz = 0(|exp-«/A|) 

(here again O (•) regards the situation when !R(z) tends to — oo with a bounded imaginary 
part). 

Notice that if A < is irrational and G Rea is finitely supported then 

/ G Res o£dz = O (| exp -*/a|) . 
J"fo(.z) 

1.2. Time spent near the singularity. In applications (for instance in the study of 
real analytical vector fields of the plane) it is important to estimate the Dulac time, that 
is the time it takes to drift from (z, w*) to (z^.Vn (z)) in the flow of the vector field. In 
the case of Xr this time is obviously 

dz 



^ J-/r(z) ^ 

Multiplying Xr by a holomorphic unit U does not change the underlying foliation (i.e. the 
Dulac map), although it does the Dulac time Tr,u- The later is obtained by integrating 
a time-form 4 . P. Mardisic, D. Marin and J. Villadelprat derived in [MMV08] the 
asymptotic expansion of the Dulac time for (deformations of) real planar vector fields 
when A is rational. We wish to complete their study in the complex setting. 
One can choose the time-form as 

^ ' ^ Xxll (x, y) ' 

so that the next result holds. 

Theorem 1.6. Take U G O (U) such that U (0, 0) ^ 0. Assume that IA is chosen in such 
a way that, in addition to conditions (X) and (R), the holomorphic function U\u never 
vanishes. Then the Dulac time is holomorphic on $7* and 

Tr,u(z) = [ f Z . 

J-jM*) XUo£ 



A meromorphic 1-form r is a time-form for a vector field X when r (X) = 1. 



COMPUTING DULAC'S MAPS OF (ALMOST EVERY) NON-DEGENERATE SINGULARITIES 6 



If moreover 5ft (A) < 0, the condition (a) holds and U — Uq £ x a yC {x, y} then, as 5ft (z) 
tends to — oo with bounded imaginary part, 

r R,u{z)~ [ , f Z p = [ +o(\zexp-*/\\) . 

J^ n (z) *U o £ J (z) A(7 Res o £ 



(z) XU ° £ J~f (z) 

The integral subtracted on the left-hand side boils down to A ^~ z ^ if U—U (0, 0) belongs 
to x a yC {x, y}. We mention that this situation can always be enforced by preparing the 
vector field: 

Lemma 1.7. Let Z be a germ of a holomorphic vector field near an isolated, non- 
degenerate singularity with ratio of eigenvalues A ^ K>o- Then there exists a £ N>o 
satisfying (a) and a choice of the local analytical coordinates such that Z = UXr for 
two germs of a function satisfying R £ x a yC {x, y} and U — U (0, 0) £ x a yC {x, y} with 
U (0,0)^0. 

This lemma is plainly trivial when A ^ M<o: in that setting Z is locally analytically 
conjugate to its linear part, corresponding to R = and U = est. When A < is 
irrational the vector field is formally linearizable and can be put in the sought form for 
any finite order a £ N, particularly one such that a + j > 0. When A = —2 for p and 
q positive co-prime integers, resonances may appear and Z may not be even formally 
orbitally linearizable. These resonances correspond to pairs (n, m) of integers belonging 
to (q,p)N (those for which nX + m = 0) in the Taylor expansion of R and U — U (0,0), 
and as such cannot appear for an index n lesser than q or for m = 0. It is thus possible 
to cancel out formally the (q — l)-jet with respect to x and the 0-jet with respect to y of 
the given functions, meaning we can take a := q. Then a + j > 0. The fact that this 
formal transform can always be chosen convergent is well known. 

1.3. Extension of the results to other singularities. Condition (X) is satisfied except 
for some cases when {A, 1 /x] n N ^ 0. The heuristic is that these singularities may not 
possess sufficiently many separatrices: the resonant node (A ^ 0) and the saddle-node 
(A = 0, exactly one nonzero eigenvalue) admit only one in general. The former case is not 
very interesting since it corresponds to vector fields which can be analytically reduced to 
polynomial vector fields (Poincare-Dulac normal forms) for which explicit computations 
are easily carried out. The geometry of the foliation itself is quite tame and completely 
understood. Save for some minor and technical complications, the framework we present 
can be adapted to encompass this case, although the trouble is not worth the induced 
lack of clarity in the exposition. 

The case of the saddle-node is richer. In [MT13] we prove that Theorem 1.1 holds in 
that case. When the saddle- node is not divergent {i.e. it admits two separatrices) it can 
be brought in the form (X) and the Dulac map admits an integral representation as in 
Theorem 1.2 (2) 

Vr (z) = (z) + (R- iix k ) o £ — 

J 7h ( z ) exp{kz) 

where (k,fi) £ N>o x C is the formal invariant of the saddle-node and T>^ x k is the Dulac 
map for the normal form which can be explicitly computed: 

_ , . . . exp (— kz) — exp (— kzA 

When the singularity is a divergent saddle-node it is possible to obtain an integral rep- 
resentation as well as a sectoral asymptotic behavior. We refer also to [LorlO] for more 
details. 
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After Seidenberg's reduction of its singularity a (germ of a) nilpotent foliation only 
has singular points either of non-degenerate or of saddle-node type. As a consequence 
the work done here and in [LorlO] is somehow sufficient to analyze more general Dulac 
maps. Yet there is a special case where it is not necessary to perform the reduction of 
singularities to be able to carry out some computations, which is in fact the most general 
formulation of the framework we introduce here, corresponding to vector fields in the 
form generalizing (X) 

X R = X + RY 

where: 

• Xq and Y are commuting, generically transverse vector fields, 

• Y admits a holomorphic first-integral u with connected fibers, 

• RE u a C {x, y} for some a > 0. 

Being given both a transverse disc £ meeting a common separatrix of X R and Y at 
some point p*, and a transverse £' corresponding to a trajectory {u = u.^} of Y, we 
can define the Dulac map of Xr joining £ to £' by lifting paths through the fibration 
(x, y) i — y u(x,y). Then, with equality as multivalued maps on £\{p*}, we have the 
implicit relation 

H o V R = H Q oV Q x exp / Rt 

Jir 

where r is some time- form of Xr and Hq a first-integral of Xq . With further work, and 
when applicable, it should be possible to derive the asymptotic behavior of T> R near p», 
as is done here. 

1.4. Structure of the article. This paper only uses elementary techniques and is con- 
sequently self-contained. Section 2 is devoted to the study of the geometry of the foliation 
and the analytic properties of the Dulac map. There are proved Theorem 1.1 and The- 
orem 1.2 (1) in the respective Section 2.2 and Section 2.3, while Theorem 1.2 (2) is 
established in Section 2.1. This paper ends with Section 3 where the explicit computa- 
tion of characteristics f 7o r z \ G o £ dz are performed for the model J-q. Yet the core of the 
section is the study of the asymptotic deviation between f Go£dz and f , ,Go£dz. 
Immediate consequences of this estimation are Theorem 1.2 (3) and the best part of The- 
orem 1.6. We end this paper with the proof of Proposition 1.5 in Section 3.2.3, completing 
Theorem 1.6. 



2. Analytic properties of the Dulac map, geometry of the foliation 

We recall that U is some polydisc centered at (0, 0) on which R is holomorphic and 
bounded, and we take (x*,y*) E U := U\ {xy = 0}. We will impose more technical 
assumptions on the size of IA in this section but they all shall be determined uniquely by 



\R\\ := sup 
u 



We will mostly use the fact that if the supremum of \R\ on U is lesser than 1 then the 
foliation has a «controlled topological type». 

Definition 2.1. All the paths 7 we use throughout the paper are C°° maps from some 
compact interval I into U. Its starting point (resp. ending point) is written 7* (resp. 7*). 



°This first-integral can be multivalued, as is the case considered above where Hq (x, y) = x 1 / x 



y- 
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(1) Let E C U\ {(0,0)} be a cross-section transversal everywhere to Tr (for short, a 
transverse to Tr) . We introduce the groupoid (E) of equivalence classes of 
paths 7 tangent to Tr with 7* G E, up to tangential homotopy (that is, homotopy 
within a given leaf of Tr) with fixed end-points. We call it the tangential 
groupoid of Tr relative to E. 

(2) The tangential groupoid of Tr relative to E is naturally endowed with a structure 
of a foliated complex surface, which can be understood as the foliated universal 
covering of Sat_7r H (E), that is the locally biholomorphic, onto map 

a : T R (E) — > Sat^, (E) 

7 7*. 

(3) The Dulac map of Tr associated with (cc*, y*) is the holomorphic function de- 
fined on 

r, : = {7er Jl (Mn{!/ = ife,i ? io}):n(7*) = i,} 

by 

v R : r, — > n- 1 (x.) 

7^7*. 

2.1. Integral representation of the Dulac map. 

Definition 2.2. A time-form of a vector field X is a meromorphic 1-form r such that 
r(X) = l. 

Because of the peculiar form of Xr one can always choose a time-form as 

da; 
Xx 

Lemma 2.3. Let E be a transverse to Tr. For given G G O (U) the integration process 



F : 7 G r fl (E) ' ^ / Gr 

ywes rise a holomorphic function whose Lie derivative Xr-F along Xr can be computed 
by considering F as a local analytic function of the end-point 7* . Then 

Xr-F = G. 

Proof. Outside the singular locus of Xr there exists a local rectifying system of co- 
ordinates: a one-to-one map ip such that iP*Xr = Jj. In these coordinates we have 
■0* (Gr) —Go tpdt. The fundamental theorem of integral calculus yields the result. □ 

Notice that X$ admits a (multivalued) first-integral with connected fibers 

H (x,y) := x~ 1/x y, 

which means that it lifts through a to a holomorphic map, still written Hq, constant 
along the leaves of Tq and whose range is in one-to-one correspondence with the space 
of leaves of Tq . 

Lemma 2.4. Let E be a transverse to Tr. The function 

H R : T R (E) — > C 

7 1 — > Hq exp / —Rt 

is a holomorphic first-integral of <j*Tr with connected fibers. 
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Proof. The fact that Hr is holomorphic on Tr (£) is clear enough from Lemma 2.3. It is 
a first integral of a* Tr if, and only if, the Lie derivative Xr ■ Hr vanishes. This quantity 
is computed as follows: 



Xr-Hr = X -Hr + R (y^j ■ H R 
Xr- f -Rt + R 

1 



£*X R = £*X + Ro£x£* [y ^ 



£*Xq — A— — h — — 
oz aw 

d\ d 



Since [y-§yy Ho = Hq our claim holds. The fact that Hr has connected fibers is a direct 
consequence of both facts that Hq also has and Hr\y: = Hq. □ 

Corollary 2.5. We have 

T>r = 2? xexp^i?r. 
Proof. For any path 7 e T* we have the relation Hr (7*) = Hr (7), that is 

H (7) cxp f -Rt = H (7*) ■ 

J-y 

The conclusion follows since 7 1— > Ho (7) is linear with respect to the y-coordinate of 7 
when x* is fixed. □ 

2.2. Foliated universal covering. Proving Theorem 1.1 amounts to proving that the 
leaves of the foliation T := £*Tr are simply-connected. Write 

d 



where 



\ dy J dw 

The vector field £* Xr is holomorphic and regular on the complex rectangle 

U := {3*0) < lnp, U(w) < lnr} . 
It induces a foliation T transversal to the fibers of 

II : (z,w) 1 — ► z, 

hence the leaf C Po is everywhere locally the graph {w = / (z)} of some unique germ of a 
holomorphic function defined in a neighborhood V (po) of IT (po). Because of this property 
the boundary of C Po is included in the boundary of the domain of study 

dU = (z) = In p or S> (w) = In r} . 

The whole argument to come relies on the existence of a family of curves included in 
C Po which project by IT on line segments of constant direction. A cycle 7 within C Po 
will therefore be pushed along those curves, as if repelled by the beam of a searchlight, 
resulting in a homotopy in C Pa with a path 7 whose projection bounds a region of empty 
interior. The compact image of IT07 can be covered by finitely many domains (V (Pn)) n <^, 
showing finally that 7 (and 7) is homotopic in C Po to a point. 
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Definition 2.6. (See Figure 2.1.) For v G {5ft (z) < hip}, < 8 < n and i) e S 1 the 
domain 

S(v,&,5) := {z : 5ft (z) < lnp, |arg(z - v) - argi? < 8} 

is called a searchlight beam of aperture 28, direction i? and vertex v. If v = II (po) we 
say it is a stability beam when the real part of the lift in C Po , starting from po, of an 
outgoing ray t > i-> v + t9, with |arg e /i?| < 8, is decreasing. 

Remark 2.7. This particularly means that the outgoing ray lifts completely in C Po as long 
as it does not cross {5ft (z) = lnp}. 




i) 



The theorem is now a consequence of the next lemma: 

Lemma 2.8. There exists 8 G]0,7r] and $ G S 1 smc/i i/iai /or a// po = (zoi^o) S W i/ie 
searchlight beam S (zo,$,S) is a stability beam. We can take 

X_ 

w 

8 G ]0,arccos(||i?||p a ) [, 

so that one can take 8 as close to f as one wishes by sujficiently diminishing p. Besides 
for any integral curve of £*Xr of the form t i— > (zo + £0, w (i)) wi/i i > 0, f 6 S 1 and 
w (0) = wq we have the estimate for 5ft (9) =/= 

exp (a5ft (zo)) 



(2.1) 



Ui (t) — Wq — t 



< 



and taking the limit 5ft (9) — > 



(2.2) 



w (t) — wait 



|A5ft(6»)|a 

exp (a5R (zq)) 



|1 -exp(a,m{6))\ 



< 



IAI 



ill-Rl 



Proof. The lift in J 7 of a germ of a ray z (i) = zo + with # G S 1 and i > 0, starting 
from po is obtained as the solution to 

w ,. s l + Ro£(z(t),w(t)) 



(*) 



A 



(0) = w , 



that is 
(2.3) 



w(t) = -(l + RoS(z + 9t,w(t))) 
A 
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The function t n- (p (t) := 5R (w (t)) is therefore solution to the differential equation 



(2.4) ip = 3? 

which particularly means that 
<p(t)-StU I 



(1 + R o £ o (z, w)) 



exp(a5ft(z )) p a 
|^j I exp {am {9)) < — \\R\ 



Exploiting the cruder estimate by taking 9 G ft exp (i [—5, 5]) we derive 

<p(t) < |^(cos<5-p Q ||i?||) <0. 

Since ip (0) < lnr it follows that if (t) < lnr as long as 5ft ((z (t))) < In p, which is our first 
claim. 

Integrating both sides of the estimate yields 



U[w(t)-w -t~ 



< 



exp (a5R (zq)) 
\\$l(9)\a 



\R\\ |1 - exp (aiSft (6»)) 



The study we just performed can be carried out in just the same way for the imaginary 
part of w, yielding 



3 ( w (t) — w — t 



< 



exp (a5R (^o)) 
|A5ft(6»)|a 



\R\\ |1 - exp (atSft (6»))| 



proving the sought estimate. 

An immediate consequence of the hicompressibihty of the foliation is the: 



□ 



Corollary 2.9. For any p G hi the leaf C po of J- passing through po is the graph of a 
holomorphic function on the simply- connected domain 



n (c 



Remark 2.10. As a closing remark we should stress that the «roughness» of dVl po is 
controlled by the aperture 28 of the stability beam, which can be taken as close to 7r 
as one wishes, and by the direction ft (which is that of the model). This is a kind of 
«conic-convexity» which forbids 9f2 Po to be too wild. In fact the closer 5ft (z) is to -co 
the closer fi po is to {5ft (w + Z ~ ZQ ) < lnr , 5ft (z) < lnp} near z. 



W(z) = hip} 
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2.3. The domain of definition of the Dulac map : proof of Theorem 1.2 (1). 

We continue to write T>r for the Dulac map of associated to some couple 
expressed in the coordinates £ (i.e. understood as a holomorphic function of (z,w)). We 
recall that for any po € U the leaf C Po of J- passing through po projects on 



We write the connected component containing z* . 

(2) // 5ft (A) > then £1 = f2» is a simply- connected domain. Particularly adh(£l») fl 
{5ft (z) = hip} is a nonempty line segment. For every N £ N>o there exists r > 
r' > such that this line segment contains at least lnp + iCy(z*) + [— 7riiV, 7riiV] 
for every 5ft (u>*) < lnr'. 

(3) 7/ 5ft (A) < there exists < p' < p depending only on a, A and \ \R\\ such that for 
every 5ft (z„) < In p' and N S N>o the domain contains some infinite half-band 
{5ft (z) < k' , | Of (z — Z*)| < irN} with k' < 5ft (z*) depending only on N, a, A and 



We first mention that is clearly open since if one can link a point (z, w* ) to (z» , 2?^ (z)) 
with a compact tangent path 7, whose image is included in the open set U, then surely 
this is again the case for a neighborhood of z. The rest of the section is devoted to proving 
the remaining items. In doing so we build an explicit tangent path linking (zo,w*) to 
(z*, T>r (zo)), see Proposition 2.12 below, which will serve in the next section to establish 
the asymptotic expansion of the Dulac map through the integral formula of Corollary 2.5. 
We underline right now the fact that the projection 7 of that path through II does not 
depend on w*, but only on zq, a, A, p, ||i?|| and z*. 

2.3.1. The integration path. We write po := (zo, w*). If 5ft (A) > and zo G £1 then both 
stability beams S(zq,'&, 8) and S (z*, 1?, 8) are included in Q Po and their intersection W 
is non-empty. Therefore zo can be linked to z* in £l po by following first a ray segment of 
S (zo, 1?, d) from zq to some point z\ in W, then from this point backwards z* along a ray 
segment of S (z*, 1?, <5), as illustrated in Figure 2.2 below. 




'/jo 



1? 




FIGURE 2.2. The path of integration 7 when 5ft (A) > 0. 
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On the contrary if 3ft (A) < the candidate region W could be beyond {3ft (z) = In p}. 
The construction must therefore be adapted. 

Proposition 2.12. Assume D?(A) < 0. There exists k£I depending only on a, X, \\R\\ 
and 2* for which the following property holds: for every zq £ one can choose a path 
7 : zq — > 2* with image inside fl po in such a way that 7 is a polygonal line of ordered 
vertexes (z Q , z±, z 2 , 2 3 , 2*) with (we refer also to Figure 2.3 below) 

• z\ = max {k, 3ft (zq)} + i3 (zq), 

• arg (z 2 - zi) = arg?9 ± S, 

• $t(z 2 ) = 3ft(z 3 ) <lnp and \% (z 3 - z 2 )\ < 2ttN + tan(|argtf| + 6) (lnp - k), 

• arg (2* — z 3 ) = arg 1? ± S . 




FIGURE 2.3. The path of integration 7 when 3ft (A) < 0. 

Remark 2.13. We could have made a similar construction without the use of n (i.e. by 
joining directly zq to some z 2 ), but we need it in order to obtain uniform bounds with 
respect to 5 (20 — z*) in the next section. 

Proof. Write 

ti± := i?exp(±i<5) . 

Either 3ft ($ + ) if 9(A) < or 3ft (■$_) if 3 (A) > is positive, let us assume for the sake 
of example that (A) < 0, the other case being similar in every respect. There exists 
K < 3ft (2*) such that arg?? = arccos (||i?|| cxp (ok)): in other words the ray segment 
(zq + R>q) n {3ft (z) < k} is included in fl po . Obviously k depends only on a, A, ||i?|| and 
2*. We take for 7 the polygonal line of ordered vertexes (20, zi, z 2 , z 3 , 2*) built in the 
following fashion. 

• If zo s {3ft (z) < k} then the partial ray (z$ + ■doM.yo) H {3ft (z) < k}, included in 
fl Po according to Lemma 2.8, leaves the region at some point z\ with 3ft (z±) = re. 
Otherwise we set z\ := zq. 

• Both rays {zi, z*} +i9-(_M>o, included in n po . intersect the line {3ft(z) = In p — e} 
in, respectively, z 2 and z 3 for e > very small. The worst case scenario to 
bound 1 3 (23 — 22)! happens when z\ = re + iS (nN + 2„), which gives the claimed 
estimate. 

The line segment [2 2 ,2 3 ], and therefore the whole image of 7, is included in Q po thanks 
to the next lemma: 
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Lemma 2.14. If 5ft (A) < then adh(£7 po ) n {5ft (2) = hip} is a nonempty line segment. 
If A ^ R>o there exists p > such that the same property holds. 

Proof. In the case 5ft (A) > we have max {5ft ($exp (±i<5))} > 0; say, for the sake of 
example, that := 7?exp (18) has positive real part. If A is not a positive number this 
property can be secured by decreasing p and taking i5 as close to ^ as need be. Take 
a path r connecting two points of adh (ft pa ) PI {5ft (2) ~ lnp} (which is a non-empty set) 
within fl po . Let / be the line segment of {5ft (2) — In p} joining those points. The ray 
p — i? + ]R>o emitted from some p £ I separates {5ft (z) < In p} into two connected regions. 
Since T starts from one of them and lands in the other one, the curve must cross the 
ray at some point q G ft Po - The ray q + t9 + R>o is included in ft po since it lies within a 
stability beam, while it contains p in its adherence. □ 

□ 

2.3.2. The dual searchlight's sweep. 

Lemma 2.15. When 5ft (A) > the searchlight's beam S (z , 8) is included in ft for 
any Zq G ft. When 5ft (A) < the beam S (z*, — 8) is included in ft. 

For any z£S (zq, — "d, 8) we can link (2, w*) to some point (zq, w) with 5ft (w) < 5ft (w*) 
by lifting in J- the line segment [z, zq\. Therefore the lemma is trivial in the case where 
5ft (A) < 0. On the contrary when 5ft (A) > the lemma is a consequence of the next one. 

Lemma 2.16. Assume that 5ft (A) > 0, zq £ Q and let 7/ := 5ft (u;»). Then for any other 
choice of w* with real part lesser or equal to 77 we have zq 6 Q as well. 

Proof. We set up a connectedness argument. Let B := {w* : 5ft (ui*) < 77} and A := 
{u>* : w t <E B and zq 6 £1}. By assumption A is not empty, and it is open in B for 
the same reason that O is open. More precisely any u>* 6 A admits a neighborhood 
V in B such that the image of 7 is included in fi( ZOiU ,) for every w G V. Let now a 
sequence (w n ) neN C A converge towards Woo £ B. If 5ft (A) > then the image of 7 is 
included in the union of the two stability beams S := S (zq, 8) U S (z*, 7?, 8) which are 
themselves included in every O( zo iWn y Because the real analytic curves defining dQ/^^-s 
vary continuously when w does we have S C £l( Zo . Wao ) also. In particular z* £ Slj^^^) 
and z belongs to £!* for := w^,. The former property implies in turn that A is a 
closed subset of B and as such spans the whole region B. □ 

Remark 2.17. When 5ft (A) < the above argument does not work since the image of 7 
must sometimes leave the (adherence of the) union of stability beams emitted by Z\ and 
2* (when it visits [22,23]). Nothing guarantees that the limiting domain ri( ZOiU)oo ) does 
not interrupt [22,23] at some point. 

2.3.3. Proof of (2). The fact that ft is simply-connected is a consequence of Lemma 2.15, 
in the way that what was previously done to prove the simple-connectedness of O po . 
From this lemma also follows the fact that adh (O) n {5ft (2) = \np} is a nonempty line 
segment, as can be seen by adapting in a straightforward way the proof of Lemma 2.14. 
In particular ft is connected. 

To prove that adh(fi) n {5ft (2) = Inp} can be arbitrarily big provided that 
be sufficiently small it is sufficient to invoke the fact that {y = 0} is the adherence of 
a separatrix of J-r, so that Tr (£) contains elements winding more and more around 
{x = 0}. 

2.3.4. Proof of (3). Because S (2*, — t?,<5) C f2* we only need to ensure that 3(7?+) and 
3(t?_) have opposite signs. This can be enforced by taking 5ft (2*) negative enough, i.e. 
by taking 8 as close to ^ as need be. 
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3. ASYMPTOTICS OF THE DULAC MAP 



Since T>r is naturally denned on the universal covering of U\ {xy = 0} we keep on 
working in logarithmic coordinates 

(x,y) = £ (z, w) = (expz,expu>) . 

We fix once and for all a preimage (z*,w*) G E^ 1 (x*,y*). For the sake of concision 
we make the convention that an object X hatted with a tilde stands for its pull-back in 
logarithmic coordinates X := £* X . Notice that the time form r is transformed into 

r = — dz . 
A 

To obtain the image of x by the Dulac map T>r we need to compute the integral 

Rt 

7rO) 

where jr is a path tangent to J-r linking (x, y*) to some point of II -1 (a;*), and we intend 
more precisely to compare this value with that of 



Rt 



7o0) 



which can be explicitly computed (Sub-section 3.2). The section's main result is the 

Theorem 3.1. Let N G N>o be given. There exists a constant M > depending only on 
N, A, a, p, \\R\\ and 8 such that for any bounded G G O (U) l~la; a C {x, y} and all x = exp z 
with 3 (z — z„) < ttN one has 



\y*\ \x\ a ■ 



Remark 3.2. Since (a + > we have \x\ a = o (jx logx|) , proving the first part 
of Theorem 1.2 (3). 

3.1. Approximation to the formal model : proof of Theorem 3.1 (1). We make 
here the hypothesis that 5ft (A) > 0. We need to compare this integral and the one obtained 
for the model, i.e. bound 







dG 


f Gt- f Gt 


< M 






dy 



A(z„) 



G (z, w R (z, zo)) - G (z, w (z, z )) dz 



7(^0) 



where 7 (zq) is a path linking zq to 2* within £l(z a ,iv*) ano - 2 ^ W R ( 2 i z o) is hs lift in Tr 
starting from (zq,u>*). We mention that 



W {z, Zq) -- 

For any (z, Wj) G U we have the estimate 



z- z 



G (z, W2) — G (z, wi) < |exp(az)| 



dG 



dy 



|expw2 — expwi 



so that 

|A(zo)| < 
Setting 



dG 



dy 



|exp (az + w (z, z )) (exp (w R (z, Zq) - w (z, z )) - 1) dz\ . 
Dr{z,z ) := \w R (z,z ) -w (z,z )\ 
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and taking exp z — 1 < \z \ exp \z\ into account we derive 



|A(z )| < 



dG 




dy 


J ^ 



/ cxpK(az + w (z, z )) D R (z,z )expD R {z,z ) \dz\ 



The proof is done when the next lemma is established: 

Lemma 3.3. There exists a constant K > 0, depending only on N, A, a, p, \\R\\, z* anc 
8, such that 



sup D R (7 (t) , z Q ) < K 
t 



where 7 is the integration path built in Proposition 2.12. The values of K is explicitly, if 
crudely, determined in the proof to come, and can surely be sharpened. 

Proof. Invoking the estimate (2.1) from Lemma 2.8 and setting 



C 2 
Ca 



a|A| 
Ci 

a P a d 



we know that, using the number k obtained in Proposition 2.12, 



sup D R (z,z ) < Ki := Ci (p a + exp (an)) 

sup Dr (z, z ) < K 2 := Ki + C 2 (p a + cxp (a«)) 

z£[zi,z 2 ] 

sup D R (z,z ) < K 3 := K 2 + C 3 (2nN + t&n{\&rgd + \)(\np~ n)) 

[22,23] 

sup D R (z,z ) < K := K 3 + C 2 {p a + cxp 3ft (az*)) . 



□ 



We conclude now the proof starting from 



|A(z )| < K 



dG 
dy 



exp 5ft (K + - z a/\) / cxp 5ft 

J ■y 



+ — I z J \dz\ 
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Let a := |a+ j\. We bound each partial integral := Jj 2 z , exp 5ft ((a + j) z) |dz| 

in the following manner: 



lo^i < exp 5ft a 



~0 



k-SR(zo) / / 1 
exp ( i5ft ( a + — 



dt 



exp 5ft ((a + 4) k) — exp 3? ((a + i) z ) 



5R(a+±) 



< 



< exp 5ft ( ( a + — ) z\ 



In p — k 



exp ( t5ft ( [ a + - ) tf H 



df 



< exp(a|2;i| 



exp (a (In p — k)) — 1 



< 



exp 



(a (V« 2 +7r 2 A^ 2 + In p - k\ ) , 



h->3 < cxp(a|z 2 



3(23-22) 



< exp [a \z 2 \ + 



exp t 



|3 (z 3 - z 2 ) 



Idtl 



< exp (a(|z*| + ttTV + tan(|arg$| + (5) (hip- k))) 
1 



x exp 
< exp I a 



(27rAT + tan(|argi9|+<5)(lnp-K)) ) , 
^ 5R (z*) 2 + 7r 2 iV 2 + lnp - 5ft (z»)^ J . 



In particular the dominant integral in the above list is io->ij so that there exists a constant 
M, satisfying the required dependency properties, with 

dG 



\A(z )\ < M 



dy 



exp 5ft (w* + qzq) . 



Since 5ft (a + j) > we have 



A(z ) = o (jz |exp5ft(^ 
as expected. 

3.2. Study of the model. 

3.2.1. Explicit computation. We want to compute for n, m G N the functions defined by 



Tn.irt (Zj 



tO) 
exp 



exp (nu + mwo (u, zq)) du 



exp ( ( n + — ) u ) du . 



If n + m /A = then A = — p/q, withp and q co-prime positive integers, and (n, m) = k (q,p) 
with fceff. In that case, and when k > 0. 

(3.1) T kq , kp {z) = (z* - z) exp (k (pw* + gz)) = O (|z exp 5ft (az)\) . 

The other case n + m /x ^ is not more difficult: 

exp ((n + m /\) (z* — z)) — 1 



^n,™ (2) = exp (mm, + nz) 



n + m /\ 
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One can see easily that as n + m /x tends to zero (which may happen if, and only if, 
A is a negative irrational) the function T„ jm grows in modulus. The dominant support 
introduced in Definition 1.4 allows to discriminate between two kind of growth rate. 

3.2.2. Resonant support. We show now that the resonant support consists of (quasi- 
resonant monomials only. 

Lemma 3.4. Assume that A < and a + i- > 0. 

(1) If X — —p/q < is a rational number then 

Res (a, A) = {k (q,p) : k £ N , kq > a} . 

(2) If X is a negative irrational we denote by (—P k /q k ) ke m its sequence of convergents. 
Then 



Res(a,A) = {(qk,Pk) ■ k £ N , q k > a} 



Proof. 



(1) Because we have n > a > q the relation \nX + m\ < ^ becomes 

\np — mq\ < — < 1 . 

2n 

Hence np = mq and since p and q are co-prime the conclusion follows. 

(2) This is a consequence of the well-known result in continued-fraction theory: if 

< i~ 2 /2 then (p, q) is one of the convergents 



of |A| 



£ Q>o is given such that 



□ 



3.2.3. Dominant terms: proof of Proposition 1.5. Nothing needs to be proved for Go so we 
assume that G expands into a power series G (x, y) = ^2 n>a m>0 G n ,mX n y m convergent 
on a closed polydisc of poly-radii at least (p + e, r + e). Because of the Cauchy formula, 
for all n, m 

\G n , m \ < C(p + ey n (r + e)- m 

where C := sup| K | =p+e _ | tf | =r+e \G(x,y)\. 
If A is not real then 



inf \nX + ml > a 13(A) . 

(n,ra)6N 2 \{(0,0)} 

Let z — be given with imaginary part bounded by Nit for some integer > and with 
real part lesser than 



/' 



3(A) 



ft (A) 



Ntt. 



By construction of fj, we have 



ft 



(y (Z, - Z) 



— (ft (A) ft (z, - z) + 3 (A) 3 (z* - z)) < 
1 



< ft ( - (z„ - z) 



so that we derive at last 



2 IAI r m o™+ SR ( 1 A) / 7 

|T„, m (z)| < W J,^ expftf-^ 



o|9f(A)| 
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Now 



GoSdz 



< 



^ ^ G n ^rnTn,rn (^) 
n~>a , m>0 

2\\\p*W^C (r + e)(p + e) 



x exp 5ft ( - — 



a|3(A)| e 2 
= 0(|exp-*/A|) , 

ending the proof for the non-real case. In fact this reasoning goes on holding even when 
A < as long as (n, m) belongs not to Res (a, A), since in that case 

\n + m /\\ > — 
2n 

has strictly sub-geometric inverse. 

Take now A negative real and G = GR es - If A = — v/q is a negative rational number 
then (3.1) provides what remains to be proved. Assume now that A is irrational. Because 
|exp z — 1| < \z\ when 5ft (z) < we have for (n, m) 6 Res (a, A) 

\T n , m (z)\ < \(z- z*)exp(nz)\r m . 

Therefore 



G R es ° £dz 



7o (z) 



< C £ (p + e)- n (r + 6)- m |r n , m (z)| 

(n,m)£Res(a,A) 

/ exp 3? (z) \ " / r 



< C\z,-z\ J2 

(n,m)£Res(a,A) 

,(r + e) (p + e) 



< C- 



|z* — z| cxp 5ft (fflz) 



0(|zexp-VA|) 



as expected. 
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